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ABSTRACT 

A logarithmic spiral is a prominent feature appearing in a majority of observed galaxies. This 
feature has long been associated with the traditional Hubble classification scheme, but historical 
quotes of pitch angle of spiral galaxies have been almost exclusively qualitative. We have developed 
a methodology, utilizing two-dimensional fast Fourier transformations of images of spiral galaxies, 
in order to isolate and measure the pitch angles of their spiral arms. Our technique provides a 
quantitative way to measure this morphological feature. This will allow comparison of spiral galaxy 
pitch angle to other galactic parameters and test spiral arm genesis theories. In this work, we detail 
our image processing and analysis of spiral galaxy images and discuss the robustness of our analysis 
techniques. 

Subject headings: galaxies: spiral; galaxies: structure; galaxies: fundamental parameters 



1. INTRODUCTION 



Approxi mately 60% of galaxies in the local Universe 
are spiral (Buta 1989). A considerable number of these 
spiral galaxies show Grand Design structure, where the 
spiral pattern is uniform and spans the entire disc of 
the galaxy. In these galaxies, the spiral pattern is of- 



method, based on a Two-Dimensional (2-D) Fast Fourier 
Transform (EFT) algorithrrj^ which decomposes images 
into spirals of different pitch angles and numbers of arms. 

A long-standing and quite successful theory of spiral 
structu re in galaxies is the quasi-stationary density wave 
model ( Lin fc Shu|1964 ). As gas enters this density wave, 



1998), and 



ten logarithmic in nature (Seigar & James 
so their appearance is scale independent, 'ine best geo- 
metric measure for logarithmic spirals is the pitch angle, 
and this can be measured for any galaxy in which spiral 
structure can be discerned, independently of the distance 
to the galaxy. 

It is worth noting that spiral structure has been ob- 
served for over 150 year^with no clear standard for 
quantitative measurement naving emerged, even though 
it correlates well with other important features of galax- 
ies, s uch as central sup ermassive black hole (SMBH) 
mass ( Seigar et al.||2008 ). Furthermore, spiral arm pitch 
angle could serve as a means to discriminate between ri- 
val theories for the formation of spiral structure in galac- 
tic discs. The aim of this paper is to present one such 



^ Spiral structure in galaxies was observed as early as 1845 (Lord 
Rosse's sketch of M51; ^Hcrschel^l859j . 
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These star forming re- 

gions, along with stars and gas, conglomerate together 
into spiral arms in the disc regions of spiral galaxies, 
with star forming regions found on the leading edges of 
arms, and dust (in the form of dust lanes) seen on their 
trailing edges. 

A recently formulated rival theory proposes that spi- 
ral arms are composed of identifiable groups of stars in 
highly eccentric and chaotic orbits, which originate near 
the ends of galactic bars. These orbits, though chaotic, 
keep the stars grouped in relatively narrow tubes known 
as manifolds, which are responsible for the observe d spi- 



ral structure ( Athanassoula et al. 2009b|a 2010). Ac 



cording to this theory, galaxies with stronger bar poten 



Th is code is publicly avail able for use at |http ://daf ix. 



uark . edu/-a ges/do«nloads .htrnl] and |http: //astro .hostTuaTrT 

edu/2DFFT/ 
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tials should have mo re open spiral structure. In a recent 
study of 27 galaxies, Marti'nez-Garcia (2012) found that 
« 60% of galaxies corroborate this theory and that galax- 
ies in which the spiral arms maintain a logarithmic shape 
for azimuthal ranges greater than 70° seem to corrobo- 
rate the predicted trend. 

Spiral galaxies are classified into three ma in types of 
spiral structure: grand design, flocculent (Elmegreen 
19811, and multi-armed. Grand design spirals are well- 
dehned two-armed galaxies and theoretical efforts have 
naturally focused on explaining these very striking pat- 
terns. Flocculent spirals are less regular with sporadic 
spiral arm segments. It has been proposed that the ori- 
gins of this kind of spiral are quite different from grand 
design spirals, the products of stochastic self-propagating 
star formation being acted upon by the differential rota- 
tion of the disc to create seg ments with the appearance 
of spiral arms ( Seiden fc Gerola)|1982( ). Multi-arm spi- 
rals have distinct spiral arms, not necessarily symmetri- 
cally spaced. It is likely that they formed as a result of 
galaxy harassment (frequent high spee d galaxy encoun- 
ters within clusters; Moore et al.||1996 ). Spiral arm gen- 
eration from external forces has been proposed to explain 
the genesis of multi-arm spiral structure in our own Milky 
Way (e.g., |Purcell et al.|[20Tl| ). 

In general, logarithmic spirals are goo d approximations 
of th e shape of galactic spiral arms (Seigar & James 
19981. Logarithmic spirals are defined in polar coordi- 



nates as 



roe 



(1) 



where r is the radius, 9 is the central angle, ro is the ini- 
tial radius when 9 = 0°, and -90° < < 90° is the pitch 
angle. The limits for the absolute value of the pitch angle 
are 0° and 90° , which produce a circle and a line, respec- 
tively. Pitch angle is defined as the angle between the line 
tangent to a circle and the line tangent to a logarithmic 
spiral at a specified radius. Small pitch angle absolute 
values are associated with tightly wound spirals and high 
absolute values with loosely wound spirals. The sign of 
the pitch angle indicates the chirality of the spiral, with 
positive pitch angles indicating clockwise outward wind- 
ing and negative pitch angles indicating counterclockwise 
outward winding (as seen from a given observer's posi- 
tion, i.e., above or below the galactic plane). 

In this paper we present a method for determining re- 
liable galactic spiral arm pitch angles. Given sufficient 
quality images, our software can reliably measure pitch 
angles by iterative 2-D FFT analyses. The paper is out- 
lined as follows: Sj2] describes the observations of the im- 
ages we use and our procedure to prepare those images 
for measurement through our software. ^ details the 
need for and the nature of our iterative adaptation to the 
FFT software, along with interpretation of its results. ^|4] 
discusses how we determine errors on our measured pitch 
angles. f|5] describes our image analysis and related tools 
for further image refinement and evaluation. Finally, in 
S|6]we present a discussion of our results and a few pos- 
sible applications of the code. 

2. OBSERVATIONS AND DATA ANALYSIS TECHNIQUES 

The galaxy images we use in this paper (unless men- 
tioned otherwise) come from the Carnegie-Irvine Galaxy 



Survey (CG^ |Ho et al.||201l[ ). This is a statisti- 
cally complete, magnitude-limited sample of 605 bright 
[Bt < 12.9 mag). Southern {6 < 0°) galaxies observed 
using the SITe2k CCD camera (with a pixel scale of 
0.259" pixel" "'^) on the 2.5 m du Pont telescope at the 
Las Campanas Observatory in Chile. The overall quality 
of the images is high, both in terms of resolution (median 
seeing ~ 1"), field-of-view (8.9' x 8.9'), and depth (me- 
dian limiting surface brightness ~ 27.5, 26.9, 26.4, and 
25.3 mag arcsec"^ in the B, V, R, and I bands, respec- 
tively). All CGS images have been oriented to have up 
as North and left as East (before we subsequently rotate 
images for deprojection purposes). In this paper, we use 
a subset of galaxies from the CGS sample in order to 
test our methods. For images not included in the CGS 
sample, we use the NASA/IPAC Extragalactic Database 
(NED^I^to acquire images. 

2.1. Two-Dimensional Fast Fourier Transformations of 
Galaxy Images 

A program called 2DFFT 



(iSaraiva Schroeder et al. 
'FT 



19941 accomplishes the 2-D Ft'l decomposition ot im 
ages. 2DFFT itself is an adaptation of the FOURNmu 



tine from Numerical Recipes in C (Press et al. 19891 



that deals with CCD (charge-coupled device) images 
The program is intended to analyze face-on or depro- 
jected galaxy orientations. The decompositi on is mod- 
eled on logarithrnic spir als. As pointed out by |Considere| 



& Athanassoula (|1988|, this method does not assume 
that observed spiral structures are logarithmic. It only 
decomposes the observed distributions into a superposi- 
tion of logarithmic spirals of different pitch angles and 
number of arms, which can be thought of as building 
blocks. This is analogous to the usual Fourier method 
of decomposing signals into a superposition of sinusoidal 
functi ons of different frequency. As per [Puerari et al.| 



by 



( 2000 1 , the amplitude of each Fourier component is given 



I{u,9)e-'^"'^+P''Uud9 (2) 



where u = Inr, r (radius) and 9 (central angle) are in 
polar coordinates, r^^^ is the inner radius, r^ax is the 
outer radius of the user-defined calculation annulus, and 
D is a normalization factor written as 



D = 



I{u,9)dud9. 



(3) 



/(m, 9) is the distribution of light of a given deprojected 
galaxy, in a (u, 9) plane, m represents the number of arms 
or harmonic modes, and p is the variable associated with 
the pitch angle (0), defined by 



tan(0) = 



Pn 



(4) 



with Pmax being the value of p with the highest amplitude 
for a given harmonic mode (see FigurefTl). As currently 
defined, the code calculates Equation [2]for < m < 6. 



^ http : //cgs . obs . carnegiescience . edu/ 
* http : //ned. ipac . caltech. e du/] 
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Fig. 1. — Fig. (left) - B-band (inverted color) image of NGC 5054. NGC 5054 is measured to have a position angle {PA) of 160° 
and an inclination angle (a) of 53.84°. Fig. [7p (right) - A{m,p) values for a deprojected B-band image of NGC 5054 with a measurement 
annulus defined by an inner radius of 160 pixels (41.4") and an outer radius of 508 pixels (132"). This indicates a peak in the three-armed 
spiral harmonic mode at pmax = 3.50. The equivalent single value pitch angle via Equation |4] is —40.60° {Note - subsequent sections and 
figures will revise this measurement with improved methods). 



Additionally, the code reports a phase angle ($) for the 
orientation of the spiral arm pattern, calculated as 



$ = arctan 



Im[A] 
Re[A] 



(5) 



where /to [A] and i?e[^] are the imaginary and the real 
part of A{p,'m), respectively. 

2.2. Image Preprocessing 

2.2.1. Deprojection 

An important step in measuring the pitch angle of 
a galaxy, regardless of the method, is to deproject the 
galaxy to a face-on orientation. This process assumes 
that a galaxy with the plane of its disc parallel to the 
plane of the sky will be circular. A circular galaxy 
with random inclination appears on the sky as an ellipse. 
Thus, a circular galaxy can be described by its position 
angle {PA; orientation of the semi-major axis in degrees 
East of North) and its axis ratio. In turn, the axis ratio 
can be further incorporated to characterize the angle of 
inclination (a) from the plane of the sky defined by 



a — arccos(fo/a) 



(6) 



where a is the semi-major axis and b is the semi-minor 
axis. Thus, an inclination angle of 0° and 90° describes a 
face-on and an edge-on galaxy, respectively. The position 
angle and axis ratio can be determined easily from images 
usi ng various programs, e.g. , SExtractor (Source Extrac- 
tor; iBertin fc Arnouts||1996D or the EL LIPSE routine in 
IRAI]'^|( |Tody|1986[ |Jedrzejewski|1987"l ). ELLIPSE works 



^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory (NOAO), which 



by iteratively fitting isophotes interactively to a galaxy 
image and reporting various parameters; most impor- 
tantly, position angle and ellipticity (1 — (b/a)). 

With the position angle and axis ratio well-defined, the 
galaxy can be readily deprojected. This is accomplished 
by rotating the image of the galaxy (see Figure [2k) by 
—PA (see Figure |2]d) and then stretching the x-axis by 
the a/b axis ratio (see Figure [21;). We use the IRAF rou- 
tines ROTATE and MAGNIFY to accompHsh the image 
rotation and stretching, respectively. This procedure ef- 
fectively aligns the semi-major axis of the galaxy with the 
y-axis on the image and then stretches the semi-minor 
axis to an equal length as the semi-major axis, thus cre- 
ating one unique radius for the galaxy and turning what 
was an ellipse into a circle. The process of deprojection 
is conducted in order to minimize errors in the result- 
ing mea sure ment of pitch angle. However, as discussed 
later in j |4.1[ precise deprojection is not necessary for the 
measurement of the pitch angle. Deprojection increases 
the range of galactic radii over which valid pitch angles 
may be measured, and thus decreases error in those mea- 
surements. As a result, our assumption of galaxies being 
intrinsically circular is not especially critical to the mea- 
surement of spiral arm pitch angle. In the case of highly 
inclined (i.e., nearly edge-on) galaxies, much of the spiral 
arms are hidden from sight and recovery of the intrinsic 
geometry via deprojection becomes increasingly difhcult 
with higher inclination angles. However, we have had 
success with deprojection on galaxies up to a = 77.47° 
for the case of IC 4831. Of course, the ability to ex- 



is operated by the Association of Universities for Research in As- 
tronomy (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 
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Fig. 2. — Fig. (left) - B-band (inverted color) image of NGC 1337. NGC 1337 displays a position angle of -30.97°, a semi-minor to 
semi-major (b/a) axis ratio of 0.218, and a corresponding angle of inclination from the plane of the sky of 77.41° as determined by Equation 
[6| FiP- pp (middle) - B-band (inverted color) image of NGC 1337 rotated by 30.97° {—PA) to align the semi-major axis with the y-Suxis. 
rig. \^ (right ) - B-band (inverted color) image of NGC 1337: the result of stretching the x-axis of Fig. pp by the a/b axis ratio (4.59), 
effectively deprojecting the image by circularizing the once elliptical shape of the galaxy. Subsequently, tne image is cropped and made 
square with the center of the galaxy at the center of the image. 



tract meaningful information from a deprojected image 
will strongly depend on the resolution of the image. For 
images of low to moderate resolution, it is unlikely to be 
able to meaningfully analyze galaxies with a > 60°. 

2.2.2. Image Cropping 

After deprojection, the next step is to determine the 
center of the galaxy. We assume that the center of the 
galaxy is the brightest region within the galactic nu- 
cleus using the IRAF routine IMCNTR to determine 
the brightest pixel location within a specified search re- 
gion. The apparent center of a galaxy is strongly af- 
fected by the interstellar extinction. As a result, differ- 
ent wavebands may yield slightly different center coordi- 
nates based on this routine. We have conducted a test 
of 10 randomly selected galaxies from the CGS sample 
and have determined that, on average, the positions of 
the central coordinates vary by a distance of 1.986 pixels 
(0.514") between B and I band images. This insignifi- 
cant discrep ancy is made even less important by further 
findings in §4.1[ which show that measurement of pitch 
angle does not critically depend on location of the precise 
center of the galaxy. The /MCiVTK-determined coordi- 
nates are then adopted as the center of the galaxy. The 
image is then cropped about the galaxy, with the cen- 
ter of the galaxy as the center of the cropped image (see 
Figure [2j;) , and the resulting image made to be a per- 
fect square, as required by the 2DFFT code, with odd- 
numbered pixel-sized edges in order to allow a unique 
median pixel location as the exact center of the image[^ 

2.2.3. Star Subtraction 

FFT image analysis is widely used for its mimicry of 
the human eye's ability to pick out symmetries and rep- 
etitions even in noisy or cluttered images. The 2DFFT 
code can measure pitch angles, in spite of the presences of 
many non-spiral features in a galaxy's image. An inter- 
esting feature is that though the code models the spiral 

As required by the 2DFFT code, the input image file must be 
in the form of a text file. We use the IRAF routine WTEXT to 
convert .fits images to .txt files. 



image as a superposition of spirals of different numbers of 
arms and different pitch angles, it measures the correct 
pitch angle even for harmonic modes where the number 
of spiral arms is incorrect. Thus, for a sufhciently low- 
noise image, it is not even necessary to correctly infer the 
number of spiral arms in the galaxy in order to accurately 
measure its pitch angle. 

When noise is introduced, this agreement in pitch an- 
gle measurement between the different harmonic modes 
is the most obvious casualty. One important source of 
noise is the presence of bright foreground stars, espe- 
cially when they are superimposed on the disc of the 
galaxy itself (see Figure p|. Nevertheless, the 2-D FFT 
of the harmonic mode with the correct number of arms 
(most commonly, to = 2) will still usually give a stable 
value, which seems to correspond to the correct pitch 
angle even when other harmonic modes show no reliable 
measure. In order to increase confidence in our pitch 
angle measurement, we reduce the noise by subtracting 
the foreground stars. Since the IRAF ELLIPSE function 
does not always work when analyzing star-subtracted im- 
ages (it fails to start if light from the center of a galaxy 
has been removed and it cannot locate the center), it 
is best to measure the ellipticity before performing star 
subtractiorj^ Often, this results in several of the har- 
monic modes coming into at least rough agreement with 
the one harmonic mode previously selected as the best 
single example. This suggests that foreground star con- 
tamination is a leading source of noise in the images and 
that star subtraction is a useful step. At the same time, 
it may not always be required fo r an accurate measure - 
ment. This is corroborated by Martfnez-Garci'a (2012"), 
who finds from a study of 27 galaxies, that the presence 
of foreground stars does not affect the value of pitch an- 
gle in general. 



^ We fit a Gaussian Point Spread Function to the bright stars 
and subtract them using the IRAF package DAOPHOT. 
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Fig. 3. — Fig. (left) - Deprojected B-band image of IC 4538 before star subtraction. IC 4538 is measured to have PA = 50° and 
! = 39.65°. Fig. (right) - Deprojected B-band (inverted color) image of IC 4538 after a Gaussian star subtraction was performed. 



2.3. Image Measurement 

Following pre-processing, the first step is to specify an 
inner and outer radius of the galactic disc. The inner 
radius is the location where the spiral arms begin, i.e., 
where the galactic bar or bulge terminates; the outer ra- 
dius is the location where the spiral arms cease, usually 
the outermost visible radius of the galaxy. Thus, an an- 
nulus is defined within which the 2DFFT code conducts 
its Fourier Decomposition. We then take this previously 
established procedure a step further by automating the 
code to measure many annuli so the final quoted pitch 
angle is not determined solely by one user-determined 
annulus. Our modifications and motivations behind the 
modifications to the code are detailed in the following 
section. 

3. PITCH ANGLE AS A FUNCTION OF INNER RADII 

The greatest source of human error lies in choosing 
an inner radius. Whereas it is seemingly easy for the 
user to visually identify the edge of the galaxy (i.e., the 
outer radius) , it is significantly more difhcult for the user 
to accurately specify the cessation of the bar/bulge fea- 
ture of a galaxy (i.e., the inner radius). Furthermore, 
slight error in specification of the outer radius has lit- 
tle ill effect, whereas slight error in specification of the 
inner radius may have significant effect. To illustrate 
this result, consider different values in outer radii: un- 
derestimation results in the full length of spiral arms not 
being measured; overestimation results in the sky being 
measured at the edge of the spiral arms. Since our FFT 
computations are luminosity biased, sky inclusion does 
not significantly affect the computations. On the other 
hand, consider different values in inner radii: overesti- 
mation results in the full length of spiral arms not being 
measured; underestimation results in a bright bar/bulge 
feature being measured in addition to spiral arms, this 



last case being the worst possible scenario. 

As a result of this observed sensitivity to inner radius 
selection, we run 2DFFT iteratively at different inner 
radii. This allows the user to specify an outer radius and 
calculate pitch angles for all possible inner radii within 
the defined outer radius. A scripting utility has been 
created to calculate pitch angles at all possible inner 
radii, given an outer radius, for a galaxy (see Figure [4]). 
Additionally, we have modified the memory allocation 
of the original code to allow for input image sizes up to 
2048 X 2048 pixels. 

The inner radius is a numerical artifact which should 
not affect the measurement of pitch angle. Therefore, we 
seek a harmonic mode in which we find a range of in- 
ner radius over which the measured pitch angle appears 
to be the most stable and consistent with the observed 
appearance of the galaxy. We aim, typically, for a mea- 
surement of pitch angle with an associated error of 2° 
to 4°. The resulting pitch angles can be plotted vs. in- 
ner radius in order to visually identify stable pitch angle 
regions as a function of inner radius beyond the influ- 
ence of a potential bar or oblate bulge feature. Stable 
regions are selected by several criteria; the stable region 
must be of the same sign (chirality) as the observed spi- 
ral arm windings in the image, it should be of the same 
harmonic mode as the visually observed number of spi- 
ral arms, there must not be any erratic fluctuation in 
pitch angle, and the region of stable pitch angles must 
be contiguous. In certain cases, the resultant pitch angle 
agrees in multiple harmonic modes, therefore, selection 
of the harmonic mode is not critically important. This 
allows us to focus on picking a stable pitch angle, even 
when the correct m value is ambiguous. To understand 
the code's behavior, we have conducted tests with very 
low noise images, artificially created logarithmic spirals. 
Not surprisingly, the code finds it trivial to measure the 
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Fig. 4. — Flow chart representing our iterative method for mea- 
suring pitch angle as a function of inner radii. 



pitch angle of such an image. For these synthetic spi- 
rals (see Figure [5]), the apparent stable regions are easily 
perceptible. However, in real galaxies, careful image in- 
spection and other techniques (see ^and its subsections) 
are sometimes required to pick out more visually elusive 
stable regions amidst the range of harmonic modes avail- 
able. 

To be clear, our method does not avoid having to 
inspect each image individually. This is a necessary 
and not totally undesirable requirement. Computer vi- 
sion methods are currently under develo pment to mea- 
sure g alactic spiral arm pitch angle (e.g., Davis & Hayes 
2012 1. The human eye has proven itself as the most reli- 
able tool for the geometric classification of galaxies (e.g 



Galaxy Zoo; Lintott et al. 2008|. The human operator 



is required only to inspect the image qualitatively for 
signs of gross error, not to re-perform any measurements 
quantitatively. 

An unavoidable side effect of deprojection is an oblate 
distortion of a galaxy's nuclear region. Initially, a nuclear 
galactic bulge is considered to be spheroidal. Therefore, a 
perfectly face-on galaxy (no deprojection required) with 
a nuclear bulge does not hinder the selection of an inner 
radius. This can be seen in the case of a synthetic log- 
arithmic spiral with a nuclear bulge component added 
(see Figures [5: and[5|;). Since a galaxy is deprojected 
according to the outer region of the galaxy and not the 
inner region, the nucleus can be distorted to an oblate 
spheroid. This creation of a non-spherically symmetric 
feature can negatively affect the calculated pitch angle 



in the innermost regions of a galaxy. 

The largest and likeliest source of error due to inner 
radii determination is when barred galaxies are mea- 
sured. Galactic bars are linear features and therefore 
have high pitch angles (0 ~ 90°). Inclusion of a high 
pitch angle feature into the measurement annulus of the 
2DFFT code results in a significant overall biasing of 
the resulting pitch angle towards the high side. This is 
always the case, because the highest practical limit for 
spiral arm pitch angle is significantly lower than the pitch 
angle produced by a galactic bar. The effect of galactic 
bars are illustrated in the example of a synthetic two- 
armed spiral with a bar component added (see Figures 
[5}i and[5[i). 

4. ERROR DETERMINATION 

The most obvious error is the variance about the mean 
pitch angle over the selected stable region in inner radii. 
The error is found by calculating the mean and standard 
deviation of the sample of pitch angles over the selected 
stable region. This standard deviation of pitch angle over 
the selected stable region is then weighted by the length 
of the stable region compared to the total length from 
the innermost spiral structure to the edge of the galaxy. 
Based on our observation from running synthetic loga- 
rithmic spirals through our code (see subsequent subsec- 
tions), reliable pitch angles are not measurable for inner 
radii selected beyond « 90% of the selected outer radius. 
At this point, too much information has been ignored for 
the code to accurately measure a pitch angle. 

In addition, it is important to consider the resolution 
of the 2DFFT code due to a discrete step size (see Fig- 
ure [6]). 2DFFT captures —50 < _p < 50 values in dis- 
crete steps of 0.25 for six harmonic modes (1 < m < 6). 
Therefore, only discrete values of pitch angle are pro- 
duced by the subsequent conversion of p — > 0. The step 
size of the discrete Fourier transform is the analog of 
the frequency stepsize in One-Dimensional (1-D) discrete 
Fourier Transforms, the smallest measurable frequency. 
This leads to a necessarily higher precision in the lower 
regime of pitch angle absolute values and in the higher 
order harmonic modes. The quantized error of the mean 
pitch angle due to the resolution of the code (see Figure 
u\ is added in quadrature to the previously determined 
standard deviation of the mean pitch angle to give a total 
error. The final error is therefore 



(7) 



where i?^ is the total pitch angle error, e™ is the quan- 
tized error for the dominant harmonic mode, <t is the 
standard deviation about the mean pitch angle, 13 is the 
distance (e.g., in pixels) from the innermost stable spiral 
structure (i.e., beyond the infiuence of a bulge or bar) to 
90% of the selected outer radius of the galaxy (0.9rmaa;), 
and A is the length (in the same units as used for /3) of 
the stable range of radii over which the pitch angle is 
averaged. Figure [8] serves as a good example of our error 
determination and its subsequent reduction by the use of 
star subtraction. 

Equation [7] reflects the fact that in our method, a hu- 
man researcher rather than a computer makes the final 
selection of pitch angle. That is, a balance of two main 
principles governs Equation [7] the fiuctuation across and 
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Fig. 5. — Synthetic spirals and their corresponding pitch angles as a function of inner radius below. The eight individual panels will be 
identified and discussed in the following top to bottom, left to right fashion. Fig. pp - One-armed synthetic logarithmic spiral (inverted 
color) with <f> = —20°. Fig. pp - Synthetic two-armed logarithmic spiral (invcrtea color) with a constant m = 2 pitch angle of —20°. 
Fig. pp - <p = —20° synthetic two-armed logarithmic spiral (inverted color) with a circular bulge component added. Fig. [5p - </> = —20° 
synthetic two-armed logarithmic spiral (inverted color) with a barred nuclear component added. Fig. - Pitch angle resuTts for Fig. [5p. 
AH harmonic modes display the correct pitch angle until 90% of the outer radius is reached. Fig. - Pitch angle results for Fig. ^3. 
The results for the even harmonic modes are essentially the same and accurately measure the correct pitch angle until Ri 90% of the outer 
radius is reached. Fig. [5|; - Pitch angle results for Fig. [s};. The vertical line at 143 pixels represents the radius of the circular bulge. The 
odd numbered harmonic modes have a systematically lower absolute value of pitch angle, but the even harmonic modes are unchanged by 
the addition of a circular bulge component. Fig. [5fi - Pitch angle results for Fig. [Sjl. The vertical lines at 27 and 85 pixels represent the 
semi-minor and semi-major axes of the bar, respectively. The odd numbered harmonic modes have systematically lower absolute values of 
pitch angle, just as in the case of the circular bulge component, but the innermost inner radii demonstrate intuitively high absolute values 
of pitch angle. 




Fig. 6. — Every possible value of pitch angle calculable by the 
2DFFT code. Pitch angles are determined from Pmax according to 
Equation |4] 



the length of a chosen stable region of pitch angle as a 
function of inner radius. Our process ensures that the 
error about the mean pitch angle is appropriate, based 
on the choices made by the user. For example, an er- 
ratic "stable" region or a short stable region will both be 
punished with appropriately high errors. Thus, a careful 
selection of stable region is required so as not to produce 
substantial errors. 

4.1. Inclination Angle and Galactic Center Position 

Errors 

The problem of making a poor choice of inner radius 
was addressed by altering the code so that it calculates a 
pitch angle for all possible inner radii. Other user-defined 
parameters have little impact on our results, and thus 
do not require such measures. The most important step 
is deprojection, which requires the user to measure the 
galaxy's inclination angle, presuming that the galaxy's 
disc is inherently circular. Tests with a synthetic two- 
armed spiral with pitch angle of —20° (see Figure [5]d) 
demonstrate that measurement of pitch angle is correct 
for any even number of arms and for inner radii up to « 
90% of the outer radius (see Figure Isf). For a one-armed 
synthetic spiral, all harmonic modes are in agreement 
(see Figure [S];). 

When the synthetic two-armed spiral (see Figure [5|d) 
is shrunk along one axis incrementally to simulate an m- 
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Fig. 7. — Every possible absolute value of pitch angle calculable 
by the 2DFFT code and their resulting average quantized error due 
to the discrete step resolution. A third-ordered best fit polynomial 
is fit to each harmonic mode (m) to interpolate error values at 
all points: y = ~4x IQ-^x^ + O.OOSSx^ - 0.0137x + 0.0234; = 
0.99994 for m = 1, -y = -2 X lO-^x^ + 0.0029x2 - 0.00843; + 0.0222; 



- 0.002x^ 



0.0064X ■ 



= 0.99997 for m = 2, 2/ = -4 X lO"' 
0.0214; _R2 = 0.99998 for m = 3, y = -1 x lO^^x'* + O.OOlSrc^ - 
0.00.54.-E + 0.0207; R~ = 0.99998 for m = 4, j/ = -9 x lO-^^x^ + 
0.0012a;2 - 0.0046x + 0.02; = 0.99999 for m = 5, and y = 
-7x 10-^a:3 + 0.001x2-0.0041a; + 0.0191; R^ = 0.99999 for m = 6. 



creasingly inaccurate deprojection, the results show that 
there is still a stable region of inner radii with the correct 
measure of pitch angle (see Figure^). Similarl y w ith a 
real two-armed galaxy, NGC 5247 (see Figure 10 d), it 
is of interest that an incorrect choice of inclination an- 
gle merely causes a gradual reduction in the length of the 
stable region over which the selected inner radii yield the 
correct pitch angle (see FigurejoJ)). Thus, deprojection is 
still an important step, but is unlikely to be a significant 
source of error when using the script, which calculates 
pitch angle for a wide variety of possible inner radii. 

Similarly, when choosing the center of the galaxy im- 
age, tests with a synthetic two-armed spiral (see Figure 
5p) and a real two-armed galaxy (NGC 5247, see Figure 
TOp) suggest that incremental errors in centering only 
gradually reduce the stable region without affecting the 
actual measure of pitch angle (provided the stable region 
of roughly constant pitc h an gle remains lengthy enough 
to be found, see Figure 11). Overall, these tests are a 



testament to the robustness of the 2DFFT algorithm. 

4.2. Bulges and Bars 

Our synthetic two-armed spiral was also used to study 
the effects of circular bulges and bars in galactic nuclei on 
pitch angle measurements. When a circular bulge com- 
ponent is added to the synthetic two-armed spiral (see 
Figure [sfc), the even numbered harmonic modes are un- 
affected7 whereas the odd harmonic modes are system- 
atically different with the lower harmonic modes being 



the worst (see Figure [5p). In contrast, when a bar com- 
ponent is added to the synthetic two-armed spiral (see 
Figure [5|i), the resulting value of the measured pitch an- 
gle is significantly increased at inner radii, with the cor- 
rect pitch angle value returning after the inner radius 
is beyond the extent of the bar (see Figure [5ji) . As an 
example, NGC 1365 (see Figure 12) displays a similar 
bar to the two-armed synthetic spiral with a bar added. 
These results confirm that circular bulges should not af- 
fect pitch angle, whereas the barred geometry can sig- 
nificantly bias pitch angle measurements towards higher 
values. We are therefore confident in the necessity of our 
efforts to systematically exclude barred nuclei from the 
pitch angle measurement annulus. 

4.3. Problems with Underlying Presumptions 

So far we have presumed that the pitch angle of a log- 
arithmic spiral is a meaningful quantity to measure in 
images of disc galaxies. Certainly there are very many 
disc galaxies for which logarithmic spiral patterns are the 
most obvious feature of the disc, as the human eye per- 
ceives it. Nevertheless, two important objections might 
be made concerning the measurement of pitch angles as 
a useful characteristic of galaxies. One is that the pitch 
angle may be different for the same galaxy when viewed 
at different wavelengths. The other is that the pitch an- 
gle might vary with the radius of the disc, in other words 
that the spiral is not truly logarithmic. 

4.3.1. The Effect of Wavelength, on Pitch Angle 

It is important to consider the possibility of different 
pitch angles arising in different wavebands of light and 
what physical processes that might imply. For instance, 
optical B-band images tend to trace the bright massive 
star forming regions of a galaxy and near-infrared (NIR) 
ima ges tend to trace the old stellar populations in galax- 
ies ( Seigar fc Ja^|1998[ [Eskridge et aL]|2002[ ) . The old 
stellar populat ion~traces the spiral density wave ( Seigar 



fc James]|1998[) . Furthermore, a spiral that appears ffoc- 
ient in the ii- 



cuient in the B-band may appear to have a weak grand 
d esign spiral in the nea r-infrared (Thornley 1996) 



Kendall et al. 



(2011) used a 1-D Fb"L' analysis on op- 
tical and MIK images of grand design spiral galaxies to 
measure their m = 2 pitch angles and concluded that 
a good correlation exists between galaxies being grand 
desi gn in the infrared and in the optical. |Seigar et al. 
(j2006j) demonstrates that a 1:1 relation exists between 
the B and NIR band pitch angles for a sample of 66 g alax- 
ies from a combination of the CGS ( Ho et al.]|2011 1 and 
the Ohio St ate University Bright Spiral Galaxy Survey 
(OSUBSGS ; 'Eskrid ge et al][2002l ). Alternatively, [Gros^ 
bol & Patsis (1998) propose a contrary view. They find 



a systematic trend of arms being tighter in bluer colors 
by investigating five galaxies in B, V, I, and K' filters. 
Admittedly, two of their five galaxies are tight spirals for 
which little or no change in pitch angle is observed, but 
it seems that more work with multiple filters is required. 

Using our method, we have r emeasured a su bset of 47 
of the galaxies appearing in [Seigar et al. ( 2006 they used 
an earlier version of this rnethod) and hav e also identi- 
fied a seemingly 1:1 relation (see Figure 13 and Table l]). 
Therefore, despite seemingly small-scale differences be- 
tween spiral arms in different wavelengths of the optical- 
NIR spectrum, the overall structure of the spiral arms. 
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Fig. 8. — IC 4538 B-band pitch angle as a function of inner radius for deprojected (PA = 50° & o = 39.65°) images before (left) and after 
(right) Gaussian star subtraction was performed, r / {D25/2) is plotted on the s econdary j-axis with D'?,^ (m ajor axis at the 25.0 B-mag/sq 
arcsec isophote) from the Third Reference Catalogue of Bright Galaxies (RC3; |de Vaucouleurs et al.ll99l| l. Fig, [gk (left) - a stable mean 
pitch angle of —17.52° is determined for the m = 3 harmonic mode from a minimum mner radius ot 3t) pixels (9.37') to a maximum inner 
radius of 208 pixels (53.9"), with an outer radius of 258 pixels (66.8"). This stable region of 172 pixels (44.5") occupies 67% of the galactic 
disc. EquationjT] yields = 3.17° with A = 172 pixels (44.5"), /3 = 196 pixels (50.8"), cr = 2.75°, and £3 = 0.47°. The final determination 
of pitch angle is therefore —17.52° ± 3.17°. Fig. [Sp (right) - a stable mean pitch angle of —17.98° is determined for the m = 4 harmonic 
mode from a minimum inner radius of 9 pixels (2^3") to a maximum inner radius of 235 pixels (60.9"), with an outer radius of 264 pixels 
(68.4"). This stable region of 226 pixels (58.5") occupies 86% of the galactic disc. Equation m yields = 2.61° with A = 226 pixels 
(58.5"), ,3 = 229 pixels (59.3"), a = 2.56°, and 64 = 0.35°. The final determination of pitch angle is therefore -17.98° ± 2.61°. This result 
is barely different from the result without star subtraction; the main difference is the redetermination of the dominant harmonic mode. 
The percent difference in mean pitch angle is 2.59% with a 17.67% reduction in error from the original. 



and thus the proposed density wave, is consistent across 
the optical-NIR spectrum. This is in opposition to the 
prediction of the density wave theory that different pitch 
angles are expected fo r spirals when observed in different 
bands ( Hozumi|2003 1. Although, the expected difference 
in pitch angle across wavelength is probably small enough 
that an extremely high precision fit would be necessary 
to falsify this prediction of density wave theory. 

From our experience, we have become accustomed to 
preferring B-band images in general due to their charac- 
teristic clarity of galactic stellar components. However, 
our comparison of pitch angles in different wavebands has 
convinced us that we can typically measure pitch angle 
across a wide range of electromagnetic wavelengths. In 
that regard, we have successfully measured pitch angles 
of galaxies in the extreme cases of far-ultraviolet and 21 
cm radio wavelength images when no other imaging data 
was available. 

4.3.2. Variable Pitch Angle with Galactic Radius 

Occasionally, spiral arms may appear to change pitch 
angle in the outer region of the disc, sometimes discontin- 
uously. These are more the exception than the rule and 
we have generally preferred to measure the inner part of 
the disc in su ch c ases or use more elaborate processing 
methods (see ^ 5.1 1 to mitigate the severity of pitch angle 



variability. Considering the case where a dichotomy ex- 



ists between the pitch angles measured in the inner and 
outer regions of a galactic disc, the code can be made to 
run iteratively for two separate regions of the galaxy and 
average the results to yield an average pitch angle for the 
disc. However, if the pitch angle results are subsequently 
used for building relationships to pr ocesses in the nucleus 
of a galaxy (e.g., Seigar et al.||2008 ), pitch angles for the 
innermost portion of a galaxy perhaps make the most 
physical sense and are furthermore not as susceptible to 
extragalactic interaction. It is also likely that the entire 
extent of a galaxy might not display logarithmic spirals. 
If so, our stable regions are selected to only highlight 
clearly logarithmic sections of spiral arms. 

To illustrate the case of measuring pitch angles of in- 
teracting galaxies, we have selected perhaps the most fa- 
mous case of interacting galaxies, M51 (see Figure 14 1). 
M51 consists of M51a (NGC 5194) and its companion 
dwarf galaxy M51b (NGC 5195). Due to M51a's clear 
interaction with M51b, its well-defined spiral structure 
is seen to depart from regularity close to the companion. 
For this case, our typical method of measuring pitch an- 
gle across the entirety of the galactic disc knowingly sam- 
ples the outer 40% of the galaxy, which is clearly seen to 
be disrupted (see Figure 14 o) . Just as in the case of iter- 
ative determination of pitch angle as a function of inner 
radius to omit interior regions, by alternatively selecting 



10 



Davis et al. (2012) 



0) 
(D 



c 

< 

x: 
o 



II . 



■50 



"If 



Face-On 

5° Inclination Anj^le Error 

20 
30 
40 

50 



Inclination Angle Error 
Inclination Angle Error 
Inclination Angle Error 
Inclination Angle Error 
Inrlination Angle Error 



_L 



100 200 300 

Inner Radius (pixels) 



400 



0.2 



r/(D35/2) 
0.4 0.6 



0.8 



50 



0) 

M 
0) 
XI 



c 
< 

O 



II 



-50 



Origina 
SExtractor 
IRAF 




200 400 
Inner Radius (pixels) 



600 



Fig. 9. — Inclination angle tests on pitch angle output for a synthetic two-armed logarithmic spiral (left) and a real two-armed galaxy 
(right). Fig. [sji (left) - Pitch angle results for different errors in inclination angle for the synthetic two-armed spiral in Fig. [sb. Even 
at a high degree of simulated inclination angle error, the mean pitch angles remain approximately the same despite a gra dually shrinking 
stable region across inner radii. Fig. [Sp (right) - Real two-armed galaxy inclination test using NGC 5247 (see Fig. llOp). Three angles 
of inclination are tested; Original - trie galaxy before deprojection, SExtractor - incorporates the deprojection according to SExtractor 
(PA = 30.4°; a = 36.97°), and IRAF - incorporates the deprojection according to IRAF (PA = 38.71°; a = 25.18°). Pitch angle results 
over their respective stable regions are within each other's error bars. I-band pitch angle results are - Original: —37.6° it 5.69°, SExtractor: 
-35.62° ± 10.96°, and IRAF: -36.16° ± 9.25°. 






Fig. 10.— Fig. pOk (left) - Deprojected {PA = 160° & a = 53.84°) B-band (inverted color) image of NGC 5054. Fig. \Idp (right) 
Deprojected {PA ^0° k, a = 28.36°) B-band (inverted color) image of NGC 5247. 
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Fig. 11. — Center error tests on pitch angle for a synthetic two-armed logarithmic spiral (left) and a real two-armed galaxy (right). Errors 
of 10, 20, 30, 40, and 50 pixels from the previously determined center are used for both. Fig. (left) - Pitch angle results for different 

errors in center determination for the synthetic two-armed spiral in Fig . Isb. As the error increases, the stable region gradually decreases, 
yet the approximate me an p itch angle remains about the same. Fig. \11\ > (right) - Real two-armed galaxy center test using a B-band image 
of NGC 5247 (see Fig. IlOp) after deprojection {PA = 38.71° & a = 25.18°) was performed. The same case is true for the real galaxy 
image; the mean pitch angle remains constant despite a decreasing stable region with increasing error. 
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Fig. 12.— Fig. \mi (left) - Star-subtracted and deprojected B-band (inverted color) image of NGC 1365: PA = 32° and a = 56.63°. 
Fig. \1S\ > (right) - Fitch angle as a function of inner radius for NGC 1365. A stable mean pitch angle of —34.81° is determined for the 
m = 2 harmonic mode from a minimum inner radius of 451 pixels (117") to a maximum inner radius of 812 pixels (210"), with an outer 
radius of 938 pixels (243"). This stretch of 361 pixels (93.5") occupies 38% of the galactic disc. Equation [t] yields = 2.80° with A = 377 
pixels (97.6"), /3 = 462 pixels (120"), cr = 1.17°, and t2 = 2.40°. The final determination of pitch angle is therefore -34.81° ± 2.80°. This 
galaxy demonstrates a large bar, approximately 34% of its outer radius. The absolute value of the pitch angle can be seen to gradually 
decrease from ip = —82.87° at an inner radius of 105 pixels (27.2") until (p = —36.03° at an inner radius of 435 pixels (113"), signaling the 
end of the bar. 
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TABLE 1 

Pitch Angle/Wavelength Comparison 



Galaxy Name 


Morphology 


a (deg.) 




B-Band 




I-Band 


m 


4> (deg.) 


m 


(deg.) 


ESO 121-026 


SB(rs)bc 


50.95 


2 


10.06 lb 1.30 


3 


11.24 it 3.05 


EbU 582-012 


SAB(rs)c 


50.21 


2 


21.45 ± 2.32 


2 


23.54 ± 3.54 


IC 2522 


SB(s)c pec 


44.77 


3 


— 26.73 it 4.69 


5 


—31.70 it 5.15 


iC 2537 


SAB(rs)c 


49.46 


4 


27.37 ± 3.84 


4 


34.11 ± 2.51 


IC 3253 


SA(s)c: 


67.05 


4 


— 17.53 ± 0.93 


4 


— 13.50 ± 3.03 


IC 4538 


SAB(s)c: 


39.65 


4 


— 17.98 ± 2.61 


3 


"1 o A f\ 1 oo 

— 18.49 ± 5.20 


iC 4831 


(K )bA(s)ab 


77.47 


2 


— 22.93 ± 2.41 


2 


"1 r\T 1 1 A 

— 16.07 ± 1.46 


NGC 150 


SB(rs)b: 


60.00 


2 


14.29 it 4.26 


1 


23.15 it 7.21 


NGC 157 


SAB(rs)bc 


50.21 


3 


8.66 it 0.89 


3 


9.32 it 1.02 


NGC 289 


SB(rs)bc 


44.77 


5 


19.71 it 1.95 


3 


12.29 it 1.36 


NGC 578 


SAB(rs)c 


50.95 


3 


16.51 ± 1.88 


2 


ri A A 1 OO "1 T 

24.24 ± 20.11 


NGC 613 


SB(rs)bc 


40.54 


3 


21.57 it 1.77 


3 


21.50 it 2.06 


NGC 908 


SA(s)c 


64.53 


3 


15.26 it 2.61 


2 


32.12 it 5.01 


NGC 1187 


bB(r)c 


41.41 


4 


— 21.96 ± 3.61 


4 


—21.55 ± 2.54 


NGC 1232 


SAB(rs)c 


28.36 


3 


— 25.71 it 5.43 


6 


—31.51 it 6.03 


NGC 1292 


SA{s)c 


64.53 


3 


— 15.89 it 2.30 


3 


— 16.34 it 5.75 


NGC 1337 


bA(s)ca 


77.41 


2 


— 16.53 ± 2.40 


3 


— 19.58 ± 2.45 


NGC 1353 


SB(rs)b: 


65.80 


4 


"1 O O 1 O O "1 

13.68 ± 2.31 


4 


13.21 ± 1.65 


NGC 1365 


SB(s)b 


56.63 


2 


—34.81 ± 2.80 


2 


—35.94 ± 4.31 


NGC 1559 


bB^sJcd 


55.25 


2 


— 26.61 ± 9.69 


2 


Ol OT 1 O O T 

— 23.87 ± 2.87 


NGC 1566 


SAB(s)bc 


36.87 


2 


— 17.81 ± 3.67 


2 


—35.73 ± 5.10 


NGC 1792 


SA{rs)bc 


60.00 


3 


—20.86 ± 3.79 


3 


oo A 1 o oo 

—20.45 ± 3.20 


NGC 1964 


SAB(s)b 


67.67 


2 


— 12.86 ± 3.49 


2 


—7.85 ± 14.61 




SB(r)b 


19.95 


3 


23.05 ± 7.90 


3 


16.91 ± 3.31 


NGC 2090 


SA{rs)c 


60.00 


4 


4.91 ± 0.56 


4 


5.19 ± 0.60 


NGC 2280 


SA(s)cd 


60.66 


4 


O 1 J T 1 O O T 

21.47 ± 2.87 


2 


13.98 ± 2.02 


NGC 2835 


SB(rs)c 


47.93 


3 


— 2o.97 ± 2.22 


3 


OT 1 T 1 O O 

— 27.17 ± 2.68 


NGC 2935 


(K )bAB(s)b 


38.74 


2 


"1 O /I 1 A TO 

— 15.24 ± 4.72 


2 


— 13.82 ± 5.32 


NGC 3052 


SAB(r)c: 


49.46 


3 


-18.45 it 1.59 


2 


—21.41 it 3.46 


NGC 3054 


SAB(r)b 


52.41 


3 


12.80 it 1.77 


3 


10.56 it 1.55 


NGC 3223 


SA(s)b 


52.41 


4 


"1 f\ OO 1 O 1 T 

— 10.92 ± 2.17 


6 


OT Tn 1 A Tn 

—27.79 ± 4.79 


NGC 3261 


SB(rs)b 


40.54 


6 


15.38 it 0.71 


6 


15.09 it 0.78 


NGC 3318 


SAB(rs)b 


57.32 


3 


35.58 it 5.53 


3 


17.85 it 4.89 


NGC 3450 


SB(r)b 


28.36 


6 


— 13.55 ± 0.31 


2 


"1 o r' T 1 Ad 

—18.57 ± 5.48 


NGC 3513 


SB(rs)c 


37.81 


1 


5.84 it 1.46 


1 


6.26 it 1.57 


NGC 3887 


SB(r)bc 


40.54 


4 


— 29.16 ± 4.82 


4 


—23.40 ± 4.08 


JNGC 4027 


SB(s)dm 


41.41 


1 


— 12.0d ± 5.47 


1 


— 10.58 ± 3.01 


IN ^UOU 


^ A /'cM-ir- 
O/A 1 b ) DC 


AA 77 


Q 
O 




Q 

o 




NGC 4050 


SB(r)ab 


47.16 


1 


-6.32 it 1.90 


1 


-6.85 it 1.63 


NGC 4930 


SB(rs)b 


34.92 


3 


30.29 it 3.45 


6 


38.86 it 2.24 


NGC 4939 


SA{s)bc 


59.34 


6 


11.48 it 1.71 


6 


13.95 it 2.41 


NGC 4995 


SAB(rs)b 


50.21 


2 


13.00 it 2.88 


6 


9.27 it 0.32 


NGC 5054 


SA{s)bc 


53.84 


3 


-25.57 it 3.72 


3 


-55.33 it 9.07 


NGC 5247 


SA{s)bc 


25.18 


2 


-31.94 it 5.75 


2 


-36.16 it 9.25 


NGC 5483 


SA(s)c 


23.07 


2 


-22.98 it 4.52 


2 


-22.31 it 8.12 


NGC 5967 


SAB{rs)c: 


53.84 


3 


18.26 it 2.05 


3 


25.42 it 3.60 


NGC 6215 


SA{s)c 


30.68 


4 


-27.43 it 5.85 


4 


-26.34 it 5.95 


NGC 6221 


SB(s)c 


44.77 


6 


-27.18 it 2.14 


6 


-27.41 it 3.19 


NGC 6300 


SB(rs)b 


47.93 


4 


-16.58 it 1.52 


4 


-16.61 it 6.19 



Note. — Co l. (1) galaxy name; col. (2) morphological type from the RC3 | |de Vau^ 
[couleurs et al.|i9 91); col. (3) inclination angle; col. (4) B-band dominant harmonic mode; 
col. (5) B-band piTch angle; col. (6) I-band dominant harmonic mode; and col. (7) I-band 
pitch angle. 

Not in |Seigar et al.| | |2006| l or plotted in Figure [Ts] 



an outer radius interior to disrupted outer regions, we can 
confine our meas urem ents to only the stable portions of 
M51a (see Figure [m);) and other similar galaxies. 

4.3.3. Flocculence 

Concerning the measurement of pitch angle of different 
types of spiral galaxies, flocculent spirals provide perhaps 
the biggest challenge. From our experience with floccu- 
lent galaxies, we find that our code most often finds them 
to have high-valued harmonic modes. Their characteris- 
tics can range from fragmented arms at best to chaos 
at worst. Elmcgreen fc Elmegreen, ( ,1987, ) defined a sys- 



tem of arm classes and descriptions to categorize spirals 
into categories with varying degrees of flocculence. They 
defined 12 arm classes (classes 10 and 11 are no longer 
in use) with 12 having the most orderly spiral structure 
and 1 the least. Galaxies with arm classes 1-4 are con- 
sidered flocculent, and those with arm classes 5-12 are 
grand design. 

We have subsequently ascertained the arm classes 
(based on blue imag es from the Palomar Observat ory 
Sky Survey) given by Elmegreen & Elmegreen ( 1987) for 
the galaxies we have measured for this paper and listed 
all available arm classes in Table [21 We have also ere- 
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B-Band Pitch Angle Abs. Value (deg.) 

Fig. 13. — Pitch angles for 47 spiral galaxies in both blue (B) 
and infrared (I) wavebands. The solid black line represents a 1:1 
ratio. The dashed blue line is a best-fit line, which is consistent 
with a 1:1 ratio within errors. 



ated two plots (see Figure 15) of pitch angle absolute 
value vs. arm class (see Figure 15 1) and pitch angle er- 



ror vs. arm class (see Figure 15 rjin order to investigate 



possible dependencies on arm classes. No clear relation- 
ship can be found from either plot, thus measurement 
of flocculent spirals do not appear to be inherently less 
precise than grand design spirals. However, our method 
of measuring pitch angle is very much dependent on the 
visual inspection conducted by the user. When initially 
inspecting images of possible candidate spiral galaxies 
for subsequent pitch angle measurement, it is more nat- 
ural to be drawn to grand design spirals. This selection 
bias can be seen in selection of our sample for this pa- 
per without prior knowledge of their arm classes. Of the 
48 galaxies listed in Table [2j only 8 galaxies are classi- 
fied as being flocculent, with the remaining galaxies all 
classified as being grand design spirals. It is our prac- 
tice to only attempt pitch angle measurement on galax- 
ies that display convincing evidence of definable spiral 
structure from image inspection. From our study of the 
COS sample thus far, we have been able to convincingly 
measure pitch angles for 62% of the spiral galaxies we 
have examined; 17% of the galaxies were rejected due 
to their high angle of inclination and the remaining 21% 
were omitted due to a lack of discernible spiral structure 
(of this 21%, among those with arm classifications from 
Elmegreen & Elmegreen (19871, 60% were classified as 
fiocculent). 

5. IMAGE ANALYSIS 

The sign of the pitch angle and the number of harmonic 
modes are very important for correct image analysis. As 
the pitch angle is calculated over all possible values of in- 
ner radii for a galaxy, it is not uncommon for pitch angle 
to vary drastically in different harmonic modes. Different 
harmonic modes will have different signs of pitch angle 



and even across one harmonic mode, sign changes may 
occur. The most apparent feature to the human eye, 
for discernable spiral arms, is the chirality of the spiral 
arms. As a result, harmonic modes that favor oppos- 
ing chirality can immediately be ruled out after a quick 
visual inspection of the image. 

For galaxies with visually distinctive spiral arms, it is 
simple enough to count the number of spiral arms by eye 
and adopt that number of arms as the correct harmonic 
mode. However, in flocculent galaxies or galaxies with 
galactic arm spurs, it maybe necessary to adopt other 
methods in selecting the correct harmonic mode. Typi- 
cally, the harmonic mode with the largest region of stable 
pitch angle across inner radii is the most valuable for our 
purposes. Nonetheless, other aspects of the code can 
lend a hand in identifying the dominant harmonic mode. 
The easiest method is by plotting the amplitude of Pmax 
as a function of inner radius. This will help identify 
the harmonic mode with the strongest amplitude over 
the largest radial range of the galaxy. For NGC 5054 
(see Figure 1 16)3)1 t he m = 3 harmonic mode is dominant 
999| ) over the outer 84.8% of the galaxy's 
several harmonic modes agree 



(Block et al 
radius! In" 



many cases 

reasonably well as to the pitch angle. In addition, pro- 
ducing an image of the Inverse EFT of a harmonic mode 
can help visually identify the "correct" harmonic mode 



(see [5.2) 



5.1. Symmetrical Component Significance 

It is a likely possibility that all of the arms of a spi- 
ral galaxy, especially galaxies with multiple arms, might 
not be perfectly symmetric. This could be the result of 
tidal disruption, galaxy harassment, etc. Whatever the 
reason, slight imperfections should be handled by the ro- 
bustness of the EFT. A common trend among galaxies we 
have analyzed is that some galaxies exhibit spiral arms 
which gradually tighten, or decrease in absolute value 
of pitch an gle, towa rd t he outer regions of the galaxy 
(see Figures 10 1 and 17 1). For galaxies with drastically 



asymmetric spiral arms or arms which demonstrate vari- 
able p itch angle, we use the method of [Elmegreen et aL| 
( 1992 1 to isolate the symmetrical component of a galaxy 
(see l^'igure 18) and then we perform a pitch angle de- 



termination on the symmetrical component (see Figures 
17 D and 19 d). Symmetric parts of galaxies are illustrated 
By making images from successive rotations and subtrac- 
tions. The procedure of Elmegreen et al. ( 1992( ) is 



S^{r,0) = {m~l)F{r,0)- 



[F{r,e)-F{r,e- 



m 

where for m > 2, Sm is the image displaying the to- fold 
symmetric part of a galaxy made from the original image 
F, and the subscript T stands for truncation, meaning 
that pixels with negative intensities are set to zero. For 
the case of a two-armed spiral galaxy, the S2 image con- 
sists of all bright features in the original image that have 
equally bright features diametrically across the galaxy. 
This procedure highlights symmetric emission, such as 
spiral arm spurs, star formation regions, etc., but it in- 
troduces spurious absorption features. For example, if 
there is a dust lane in only one arm, then only the bright 
part of that arm will appear in both arms of the S2 im- 
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Fig. 14.— Fig. \14^ (left) - Deprojected {PA = -154.1° & o = 19.09°) B-band (inverted color) image of M51 (M51a and its dwarf 
companion galaxy M51b) acquired from NED (imaging from KPNO 2.1 m CFIM with a pixel scale of 0.305" pixel" '^). Fig. (middle) 
- A stable mean pitch angle of 19.13° is determined for the m = 2 harmonic mode from a minimum inner radius of 54 pixels (16.5") to 
a maximum inner radius of 229 pixels (69.8"), with an outer radius of 653 pixels (199"). This stretch of 175 pixels (53.4") occupies 27% 
of the galactic disc. Equation myields = 4.76° with A = 175 pixels (53.4"), /3 = 534 pixels (163"), ct = 1.54°, and €2 = 0.78°. The 
final determination of pitch angle is therefore 19.13° it 4.76°. Due to the interaction with its co mpanion galaxy, M51a shows a significant 
departure from a constant pitch angle in the oute r re gions of the galaxy. This is seen in Fig. |14^ and at the noticeable sign change in this 
plot at an inner radius of 389 pixels (119"). Fig. (right) - A stable mean pitch angle of 16.26° is determined for the m = 2 harmonic 
mode from a minimum inner radius of 54 pixels (To. 5") to a maximum inner radius of 276 pixels (84.2"), with an outer radius of 389 
pixels (119"). This stretch of 222 pixels (67.7") occupies 57% of the measurement annulus. Equation [t] yields = 3.20° with A = 222 
pixels (67.7"), /3 = 296 pixels (90.3"), a = 2.36°, and £2 = 0.57°. The final determination of pitch angle is therefore 16.26° ± 3.20°. This 
alternate pitch angle measurement isolates the inner portion of the galaxy out to the clear break from constant pitch angle seen in Fig |14| 3. 
As a result, the unstable outer portion of the galaxy has been ignored and a more accurate pitch angle has been determined for the purer 
inner structure of this interacting galaxy. 




5 10 5 10 

Arm Class Arm Class 

Fig. 15. — Pitc h angles (left) and their associated errors (right) sorted into their respective arm classes (classes 10 and 11 are no longer 
in use). Fig. \15\ i (left) - Pitch angles (black crosses) arranged into their arm classes wit h bin ned averages (red squares connected by red 
line segments) . No clear trend is recognizable between pitch angle and arm class. Fig. | j5f ) (right) - Pitch angle errors (black crosses) 
arranged into their arm classes with binned averages (red squares connected by red line segments). No clear trend is recognizable between 
pitch angle errors and arm class. 
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TABLE 2 
Pitch Angles and Arm Classes 



Galaxy Name 


Morphology 


m 


<P (deg.) 


Band 


bource 


Arm Class 


IC 2522 


SB(rs)bc 


3 


-26.73 ± 4.69 


B 


1 


5 


IC 2537 


SAB(rs)c 


4 


27.37 ± 3.84 


B 


1 


2 


M51a 


SA{s)bc pec 


2 


16.26 ± 2.36 


B 


3 


12 


NGC 45 


SA(s)dm 


3 


-32.13 ± 3.85 


B 


1 


1 


NGC 150 


SB(rs)b: 


2 


14.29 ± 4.26 


B 


1 


12 


NGC 157 


SAB(rs)bc 


3 


8.66 ± 0.89 


B 


1 


12 


NGC 210 


SAB(s)b 


2 


-15.81 ± 3.25 


B 


1 


6 


NGC 289 


SB{rs)bc 


5 


19.71 ± 1.95 


B 


1 


12 


NGC 578 


SAB(rs)c 


3 


16.51 ± 1.88 


B 


1 


9 


NGC 598 


RAfsVd 


2 


— 33 90 + 5 72 


6450 


2 


5 


NGC 613 


SB{rs)bc 


3 


21.57 it 1.77 


B 


1 


9 


NGC 895 


SA(s)cd 


2 


— 38.50 di 4.77 


I 


1 


9 


NGC 908 


SA(s)c 


3 


15.26 di 2.61 


B 


1 


9 


NGC 1042 


SAB(rs)cd 


4 


39.50 di 4.48 


R 


1 


9 
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NGC 2980 




4 


91 47+9 ft7 


g 


I 


9 


NGC 2442 


SAB(s)bc pec 


2 


1 4 95 + 4 20 


V 


1 


7 


NGC 2835 


SB{rs)c 


3 


— 23.97 ± 2.22 


g 


1 


9 


NGC 2935 


CRMSABI's'lh 


2 


— 15.24 ± 4.72 


g 


1 


8 


NGC 3052 


RABfrV- 


3 


— 1 8 45 + 1 59 


g 


1 


9 


NGC 3054 


RABfrib 


3 


1 2 SO + 1 77 


g 


1 


9 


NGC 3450 


RB(r)b 


g 


— 1 3 55 + 31 


g 


1 


9 


NGC 3513 


RB(rs)c 


1 


5.84 ± 1.46 


B 


1 


12 


NGC 3783 


(R')SB(r)ab 


2 


10.71 ± 0.64 


B 


1 


9 


NGC 3887 


RB{r)bc 


4 


-29.16 ±4.82 


B 


1 


2 


NGC 3938 


RA(s)c 


4 


-22.37 ± 7.21 


B 


3 


9 


NGC 4027 


RB(s)dm 


1 


-12.06 ± 5.47 


B 


1 


4 


NGC 4030 


RA(s)bc 


3 


23.48 ± 5.76 


B 


1 


9 


NGC 4050 


RB{r)ab 


1 


-6.32 ± 1.90 


B 


1 


5 


NGC 4321 


RAB(s)bc 


5 


21.81 ± 3.57 


R 


3 


12 


NGC 4939 


RA(s)bc 


6 


11.48 ± 1.71 


B 


1 


12 


NGC 4995 


RAB(rs)b 


2 


13.00 ± 2.88 


B 


1 


6 


NGC 5054 


RA(s)bc 


3 


-25.57 ± 3.72 


B 


1 


5 


NGC 5085 


RA(s)c 


2 


-11.32 ± 1.77 


4680 A" 


4 


2 


NGC 5236 


RAB(s)c 


6 


-16.04 ± 1.74 


B 


1 


9 


NGC 5247 


RA(s)bc 


2 


-31.94 ± 5.75 


B 


1 


9 


NGC 5861 


RAB(rs)c 


2 


-14.91 ± 0.83 


V 


1 


12 


NGC 6215 


RA(s)c 


4 


-27.43 ± 5.85 


B 


1 


12 


NGC 6300 


RB(rs)b 


4 


-16.58 ± 1.52 


B 


1 


6 


NGC 7793 


RA(s)d 


2 


13.91 ± 4.40 


B 


1 


2 



Note. — C"^-.. g alaxy name; col. (2) morphological type from the RC3 ( |de| 
[Vau couleurs et al.lll991|l; col. (3) dominant harmonic mode; col. (4) pitch angle; col. 
(5) waveband/ wavelen gth; col. (6) telesc ope / survey imaging source; and col. (7) arm 
class from |Elmegreen fc Elmegreen llfl987]|. Source (1) CGS; source (2) Palomar 48 inch 
Schmidt; source (3j KPiNU 2.1 m cViiVi; and source (4) UK 48 inch Schmidt. 
^ 103aE emulsion. 

^ In addition to spiral arms in the disc of the galaxy, NGC 1097 displays rare m = 2 
nuclear spiral arms in the bulge. These arms display an opposite chirality to the disc 
arms with if> = -30.60° ± 2.68°. 
IllaJ emulsion. 
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Fig. 16.— Fig. \16^ (left) - Star-subtracted and dcprojccted (PA = 160° & o = 53.84°) B-band (inverted color) image of NGC 5054 
(see Fig. ITOb), overlaid with the contours of the Inverse 2-D FFT for the m = 3 harmonic mode (in red) conducted with an inner radius 
of 160 pixels (41.4") and an outer radius of 508 pixels (132"). The contours are the real part of the complex spatial function of Equation 
[9] with m = 3. The contours illustrate the different levels of amplitude for the m = 3 harmonic mode. The Inverse 2-D FFT displays a 
smgle value pit ch a ngle of —40.60° (as shown in Fig. [iji). However, the pitch angle can be seen to tighten (decrease) as the inner radius 
increases. Fig. \16\ > (right) - Plot of the amplitude of Pmax as a function of inner radius for NGC 5054, indicating the m = 3 component 
as the dominant harmonic mode for the galaxy. The m = 3 harmonic mode is dominant from an inner radius of 77 to 456 pixels (19.9" to 
118"), constituting about 75% of the galaxy's radius. 



age; this gives the false impression that there is a dust 
lane in the other arm also. This method appears to de- 
crease our error estimates when performed. For galaxies 
with apparent initial symmetry, the mean pitch angle 
is not changed significantly; e.g., NGC 5247 (see Figure 
19), the percent difference in mean pitch angle is 10.48% 
with essentially the same error. This tool seems most 
useful for galaxies that display variable pitch angle (see 
Figure 17). Error estimates can be reduced drastically 
for these cases; e.g., NGC 5054, for which the percent 
difference in mean pitch angle is 4.19% with a 72.03% 
decrease in error. Moreover, this process can act as an 
effective substitute for star subtraction. 

5.2. Two-Dimensional Inverse Fast Fourier Transform 

One of the most powerful tools provided by 2DFFT is 
the ability to run an Inverse FFT. After having depro- 
jected the images and identified the dominant harmonic 
modes, we can calculate t he inv erse of the transforms ac- 
cording to Seigar et al. (2005). The inverse transform 



can be written as 



where 

Sm{u) 



S{u,9) 



D 



Sm{u)e'' 



(9) 



g2n4jr2 



P+ 



G,n{p)A{p,m)e'P^dp. (10) 



Gmi p) is a high-frequency filter used by Puerari & Dot- 
tori p992| ) . For the logarithmic spiral governed by Equa- 



tion|4j it has the form 
G™(p) 



= e 2^ 



(11) 



This filter is also used to smooth the A{p, m) spectra at 
the interval ends {p = — 50 and p+ = 50 with dp = 0.25) 
(Puerari & Dottori 1992| ). Equation[9]is designed as such, 
to allow the user to create an inverse transform for a se- 
lected number of harmonic components. For example, 
the inverse transform can be calculated for one compo- 
nent, e.g., m = 2, or any number of components can be 
combined to yield a composite result, e.g., m = 2, 3, & 4. 

Once an Inverse FFT is created, it can be directly com- 
pared to the deprojected image of the galaxy, allowing 
us to effectively observe what the code is seeing. Figure 
16 1 and Figure [20] show images of spiral galaxies overlaid 
with contours representing the results of Inverse FFTs of 
the same galaxy. The contours are the real part of the 
complex spatial function of Equation |9] The use of these 
images to analyze a galaxy can lead to more confident 
determination of pitch angle. 

6. DISCUSSION AND FUTURE WORK 

Our modified version of 2DFFT is a powerful tool 
for accurately measuring galactic spiral arm pitch an- 
gle. Our software, combined with careful image and data 
inspection, comparative pitch angle selection, and a self- 
regulating error determination allows for reliable pitch 
angle measurements. We hope that quantitative deter- 
minations of spiral arm pitch angle will aid in galaxy 
classification, in the indirect study of central black hole 
masses and more generally in our understanding of galac- 
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Fig. 17. — Fig. \17[ i (left) - B-band pitch angle a s a function of inner radius for NGC 5054, after star subtraction and deprojection 
{PA = 160° & a = 53.84°) were performed (see Fig. |10^ ). A stable mean pitch angle is difficult to define since the pitch angle is seen to 
continually decrease from an inner radius of about 20U to one of about 425. A rough mean pitch angle of —24.52° (a significant departure 
from the single value, non-iterative measurement of —40.60°, see Fig. ^3) is determined for the m = 3 harmonic mode from a minimum 
inner radius of 123 pixels (31.9") to a maximum inner radius of 434 pixels (112"), with an outer radius of 508 pixels (132"). This stretch 
of 311 pixels (80.5") occupies 61% of the galactic disc. This measurement has a significant error due to the unstable pitch angle. Equation 
[rjyields = 12.84° with A = 311 pixels (80.5"), 13 = 334 pixels (86.5"), cr = 11.92°, and £3 = 0.92°. The final determination of pitch 
angle is therefore —24.52° ± 12.84°. Fig. \n\ > (right) - For the pure symmetrical component of NGC 5054 (see Fig. |18^ ), the B-band pitch 
angle (same deprojection parameters as Fig. |17^ ) as a function of inner radius is far more stable. A stable mean pitch angle of —25.57° 
is determined for the m = 3 harmonic mode trom a minimum inner radius of 91 pixels (23.6") to a maximum inner radius of 253 pixels 
(65.5"), with an outer radius of 593 pixels (154"). This stretch of 162 pixels (42.0") occupies 27% of the galactic disc. Equation [t] yields 
= 3.72° with A = 162 pixels (42.0"), P = 443 pixels (115"), a = 1.31°, and t2 = 1.00°. The final determination of pitch angle is 
therefore —25.57° it 3.72°, a percent difference of 4.19% in mean pitch angle with a 72.03% reduction in error from the original. 



tic morphology and its evolution. One important advan- 
tage of this means of describing galaxies is its relative 
ease of acquisition, since only imaging data is required 
to measure it. Also, it provides us with great opportu- 
nity to t est competing theories behind galactic spiral arm 
genesis ( Martfnez-Garci'a]|2012 ) . 



6.1. Comparison to Other Methods 
Our method adds one dimension to re cently published 



FFT methods (e.g., [Kendall et aLpOll ). In doing so, we 
are able to effectively use more of the inherent informa- 
tion in the images. Alternatively, 1-D methods identify 
radial peaks in intensity azimuthally about the galactic 
center by use of radial search segments that ultimately do 
not utilize the full resolution of the image. Admittedly, 
our method comes at a higher price in terms of com- 
puting power to analyze the full resolution of an image, 
but with modern computational power, this computa- 
tional expense is trivial and is easily handled by modern 
processors. Ultimately, users of both 1-D and 2-D FFT 
methods are still obliged to visually inspect images. As 
for flocculent galaxies, 1-D FFT methods are admittedly 
only equipped to handle grand design spirals. Our 2- 
D FFT, though it may occasionally encounter trouble 
with high degrees of flocculence, will perform adequately 
with flocculent galaxies and with the support of addi- 
tional image analysis methods (see ^ and its subsec- 



tions) and multi-wavelength imaging, it can confidently 
approach any galaxy with hints of spirality. Addition- 
ally, our pitch angle measurements are unique in that fact 
that we measure and quote pitch angles for multi-armed 
galaxies (dominant harmonic modes m > 2). Other re- 
searchers may have analyzed and discussed the influences 
from higher order harmonic modes, but in the end al- 
ways publish pitch angles resulting from harmonic modes 
m < 2. 

In order to compare the results of our method to other 
published methods, we have chosen a well-studied sample 
of galaxies whose pitch angles have been determined in- 
dependently, using different techniques in the literature 
For this sample, we have selected the r esults of iMartinez- 
Garci'a (2012), Kendall eTaT] ([Ml]) jMi ( |200ip , |Gros^ 
|bol & Patsis^ (1998) , and" Kennicutt ( 1981[ ) as references 
with which to compare our measurements (see Table |3]). 
These five references provi de a n ice spread in measure- 
ment methods. Kennicutt ( 1981 ) used Ha photographic 
plates and subsequent by-hand geometric measurements 
of nearby Sa-Sc galaxies to calculate average pitch an- 
gles determined from the two main arms in each galaxjP] 



* [Savchenko fc Reshetnikov| l |201l[ ) remeasured pitch angles for 
46 ot ttie spiral galaxies measured b y [Kennicutt 1 1 |1981 ') with two 
distinct methods: by-hand and 1-D FFT analysis. Kesults of both 
methods show good mutual agreement with the mean differences 
between measurements less than a few degrees in both cases. 
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Fig. 18. — Fig. \18[i ( left) - m = 3 symmetrical component (inverted color and with the same deprojection parameters as Fig. |TOk) of 
NGC 5054 (sec Fig. |iO^ ). Fig. \l8f} fripht) - m = 2 symmetrical component (inverted color and with the same deprojection parameters as 
Fig. [l0|3) of NGC 5237 (see Fig^TOjj). 




Grosbol & Patsis (1998) utilized accurate surface pho- 
tometry and derivation of axisymmetric components to 
calculate m — 2 pitch angles for five galaxies from the 
residuals of intensive image processing and Fourier com- 
ponents of the azimuthal intensity variations; Maj( 2001 1 
visually selected points along spiral arms in CCD im- 
ages of galaxies and fit logarithmic spirals to the points; 
Kendall et al. (20111 employed 1-D FFT decomposition 
to calculate m = 2 pitch angles for a sample of grand 
design spiral galaxies; and Martmez-Garcia (20121 used 
both a "slope method" and adopted a 2-D I'll algo- 
rithm similar to our own, but without our analysis of 
pitch angle as a function of inner radius. 

Table |3] compares 38 of our pitch angle measurements 
against available measurements from five independent 
sources. The mean difference between the measure- 
ments for the same galaxies are as follows: A(f>i = 
—5.13° ± 19.41° (the differe nce between our measu red 
pitch angles and those of Martmez-Garcia 2012) or 
A01 — —0.11° ± 7.38° if the outlying measurement of 
NGC 499 5 is disregarded, Ad ^, = 1.73° ± 3.58° (the 
same for IKendall et al.|[20TT|, A03 = 0.66° ± 9.97° 



(the same with the average arm measurement fr om Ma 



20011), A04 
Patsis 199^ 
KennicuttJ li 



2.93° ± 1.71° (the same for .GrosboTX 
and A(t)5 = 5.15° ± 8.66° (the same tor 
981). These differences are comparable to 



the mean error in our measurements for this sample: 
(j)E = 3.34° ± 1.94°. 

6.1.1. NGC 7083 

For individual measurements, several differences can 
be explained by our selection of a different harmonic 
mode from that chosen by the other group. For example, 
our measurem ent of NGC 7 083 differs from the measure- 
ment by ,Grosbol fc Patsis] ( |1998) ; -19.44° ± 3.21° and 



— 15.0° ± 1.0°, respectively. This can be explained by 
our selection of the m = 3 harmonic mode and thei r se - 
lection of the m — 2 harmonic mode (see Figure 21). 
In short, we believe this is a three- armed galaxy, not 
two-armed, a finding supported by the strength of our 
code's m — 3 harmonic mode and by visual inspection. 
It is of note that we find our measurement of the m = 4 
harmonic mode's pitch angle to be —15.38° ± 2.97° (see 
Figure |21| :), which coin cides with the measurement of 
[G"rosbora"Patsis| ( |1998[ ). However, our measurement of 
the m = 2 harmonic mode is not possible due to its chiral 
instability. Despite the nice agreement of pitch angle be- 
tween our m = 4 pitch angle and their m = 2 pitch angle, 
we find strong indications that the dominant harmonic 
mode is m = 3 (see Figure 21 j). 



6.1.2. NGC 4995 
NGC 4995 depicts a significant outlie r from the pitch 



angle measurement of 



Martmez-Garcia 



(2012 



78.3 



-1-5.4 
8.3 



degrees compared to oiir measurement of ia.uO° ± 2 
Their enormously high measured pitch angle is most sen- 
sibly unphysical, along with any measurement of 4> ^ 60° . 

6.1.3. NGC 1365 

A big discrepancy can be seen in comparison of our 
measurement of NGC 1365, -34.81° ±2.80° . to th e pitch 



angle absolu te value measurem ents of Ma (20011, 13 
& 17.8°, and |Kennicutt| ( p8T| ), 18° ±3°. We find strong 
indications both from our code (see Figure 12) and vi- 
sually (see Figure 22 1 that the pitch angle oeyond the 
large bar is on the high side. In order to visually com- 
pare the fit of two logarithmic spirals with different pitch 
angles, the scaling must be adjusted. According to Equa- 
tion [T] the radius of a logarithmic spiral with a higher 
pitch angle will grow much more rapidly than a logarith- 
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Fig. 19. — Fig. [Tgip , (left) - B-band pitch angle as a function of inner radius for NGC 5247 after deprojcction {PA = 20° & a = 28.36°) 
was performed (see t ig. llOp). A stable mean pitch angle of —28.76° is determined for the m = 2 harmonic mode from a minimum inner 
radius of 1 pixel (0.259") to a maximum inner radius of 331 pixels (85.7"), with an outer radius of 565 pixels (146"). This stretch of 330 
pixels (85.5") occupies 58% of the galactic disc. EquationlTlyields = 5.73° with A = 330 pixels (85.5"), 13 = 508 pixels (132"), cr = 3.56°, 
and 62 = 1.70°. The fin al d etermination of pitch ang le is therefore -28.76° ± 5.73° . Fig. (right) - For the pure symmetrical component 
of NGC 5247 (see Fig. [l8p), a stable mean B-band pitch angle (same deprojection parameters as Fig. |19^ ) of —31.94° is determined for 
the m = 2 harmonic mode from a minimum inner radius of 9 pixels (2.33") to a maximum inner radius ol 235 pixels (60.9"), with an outer 
radius of 486 pixels (126"). This stretch of 226 pixels (58.5") occupies 47% of the galactic disc. EquationlTlyields = 5.75° with A = 226 
pixels (58.5"), /3 = 428 pixels (111"), a = 2.83°, and t2 = 2.06°. The final determination of pitch angle is therefore -31.94° ± 5.75°, a 
percent difference of 10.48% in mean pitch angle with essentially the same error as the original. As a characteristic example, NGC 5247 
displays a similar pattern of agreement among even harmonic modes as the two-armed synthetic spiral with an added symmetrical bulge 
component (see Figure [5t). 



mic spiral with a lower pitch angle. In order to allow 
the radius of the cj) = 16.5° logarithmic spiral to grow at 
the same rate as the three higher pitch angle logarithmic 
spirals in Figure |22[ we multiplied its resulting radius by 
a factor of 2.8. Wnen optimally scaled, our high pitch 
angle measurement and their low pitch angle measure- 
ments can be brought into rough agreement. The diffi- 
culty in accurately measuring pitch angle increases as the 
amount of angular wrapping around a galaxy decreases, 
i.e., Omax — f for NGC 1365 whereas the spirals do not 
become significantly different until ^ ^ at the edge of 
Figure [22] Furthermore, N GC 1365 is a well-docurnented 
case of variable pitch angle ( Ringermacher & Mead 2009 1 
with high pitch angle near the arm- bar junction and low 
pitch angle in the outermost regions of the galaxy. Our 
code correctly identifies this high pitch angle near the 
arm-bar junction and provides us with our desire d inner- 
most stable pitch angle of the galaxy (see ^ 4.3.2). 



6.1.4. NGC 3513 

NGC 3513 demonstrates a case where pal (120011) mea- 
sures dramatically different pitch angles between two 
arms (18.2° & 8.1°). Our measurement of 5.84° ± 1.46° 
for the TO = 1 component is in near agreement with the 
smaller of Ma's two measurements. From inspection of 
the image of NGC 3513, it appears to have one arm with 
near-constant pitch angle and another, more ambiguous 



arm, with a clear kink in it. The forced inclusion of the 
non-logarithmic arm might also be responsible for mak- 
ing Martfnez-Garci'a's 2-D FFT measurement of 24:.2'^q'1j 
degrees even higher than Ma's individual measurement 
of the arms. 

6.2. Pitch Angle - SMBH Relation 

Strong evidence su ggests that SMBHs reside in the nu- 
clei of most galaxies (Kormendy fc Richstone||1995 Ko- 
rmendy & Gebhardt 2001). Additionally, it has recently 
been shown that a correlation exists between the pitch 
angle of spiral arm s and SMBH mass in disc galaxies 
( Seigar et al.||2008 ). The relation is such that more mas- 
sive SMBHs reside in galaxies with low pitch angle spiral 
arms (i.e., those that are tightly wound) and the least 
massive SMBHs are found in galaxies with high pitch 
angle spiral arms. Given that a significant fraction of 
galaxies in the Univer se have spiral or barred spiral mor- 
phologies ( Buta|1989 ), we wish to improve upon existing 
methods for measuring spiral arm pitch angle in order 
to quantify their structure. This measure can in turn be 
used to estimate the central black hole mass. 

One of the current widely used relationships to SMBH 
mass is stellar velocity dispersion of the bulge /sp heroid 
( [Gebhardt et al.||2000a[ [Ferrarese fc Merritt||2000| . This 
technique requires spectroscopy ot the galactic nucleus. 
Pitch angle determination only requires optical imag- 
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Fig. 20. — Fig. [goj i (left) - Deprojected {PA = 50° & a = 39.65°) B-band inverse image of IC 4538, before star subtraction, overlaid 
with the contours ot the Inverse FFT for the m = 3 harmonic mode (in red), conducted with an inner radius of 54 pixels (14.0") and an 
outer radius of 258 pixels (66.8"), demonstrating a pitch angle of —19.98°. The contours are the real part of the complex spatial function 
of Equation |9] with m = 3. The contours illustrate the different levels of amplitude for the m = 3 harmonic mode. The overlaid In vers e 
FFT tends to track the bright foreground stars. Fig. \20p (right) - Star-subtracted and deprojected (with the same parameters as Fig. [20b) 
B-band inverse image of IC 4538, overlaid with the contours of the Inverse FFT for the m = 4 harmonic mode (in red), conducted with an 
inner radius of 54 pixels (14.0") and an outer radius of 264 pixels (68.4"), demonstrating a pitch angle of —15.42°. The contours are the 
real part of the complex spatial function of Equation |9] with m = 4. The contours illustrate the different levels of amplitude for the m = 4 
harmonic mode. In the absence of bright foreground stars, all four visible arms of the galaxy are tracked by the Inverse FFT overlay. 



r/(D,,/2) r/(D,y2) 
0,2 0.4 0.6 0,8 0,2 0-4 0,6 0,8 




64 8" Inner Radius (pixels) Inner Radius (pixels) 



Fig. 21.— Fig. \2l^ (left) - Star-subtracted and deprojected [PA = 5° & a = 53.84°) B-band inverse image NGC 7083. Fig. (middle) 
- Plot of the amphtude of pmax as a function of inner radius for NGC 7083, indicating the m = 3 component as the dominant harmonic 
mode for the galaxy. The m = 3 harmonic mode is dominant from an inner radius of 145 to 384 pixels (37.6" to 99.5"), constituting about 
61% of the galaxy's radius. Fig. \21\ : (right) - A stable mean pitch angle of —19.44° is determined for the m = 3 harmonic mode from a 
minimum inner radius of 143 pixels (37.0") to a maximum inner radius of 319 pixels (82.6"), with an outer radius of 390 pixels (101"). 
This stretch of 176 pixels (45.6") occupies 45% of the galactic disc. Equation m yields = 3.21° with A = 176 pixels (45.6"), /3 = 208 
pixels (53.9"), a = 2.67°, and £3 = 0.58°. The final determination of pitch angle is therefore -19.44° ± 3.21°. 
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Fig. 22.— NGC 1365 (see Figure \12\i) overlaid with a m = 2 
spiral with = —34.81° (solid red Tmes) representing our best- 
fit pitch angle measurement (see Figure 1 12^ 3), <f> = —32.01° (short 
dashed red lines) representing our lower limit fit, <f> = —37.61° (long 
dashed red lines) representing our upper limit fit, and cp = —16.5° 
(alte rnating sho rt-lo ng dashed blue lin es) representing the average 
fit of [Ma] pOOTI l and IKennicutt] | [T98l| . 



ing. Compared to simple optical imaging, which is 
widely and readily available, spectroscopy is time in- 
tensive. Other methods such as Reverberation Map- 



ping (Gebhardt et al. 2000b) require long-term cam 
paigns to obtain multi-epoch spectra and require sig- 
nificant telescope time and allocation. Even black hole 



estimates from single-epoch spectra fVestcrgaard 2002 1 
require spectroscop y. Other tec hniques such as bulge lu- 
minosity estimates ( Kormen dy 1993; .Kormendy fc Rich- 
stone|1995||MaKorrian et al. [1998,, Marconi fe Hunt 
Ri: 



Haring & Rix 2004p require bulge decomposition. One 



spec ific bulge lum inosity estimate incorporates Sersic In- 
dex (Sersicll963) measurements of elliptical galaxies and 



the bulges of disc galaxies and relates them to SMBH 
mass ( [Graham fc Driver ^^2007| . 

6.3. Evolution of Pitch Angle with Redshift 

No matter how a logarithmic spiral is scaled, pitch an- 
gle is unaffected. This allows pitch angle measurements 
for distant galaxies to be considered equally valid as those 
for local galaxies. For distant galaxies, as long as spiral 
arms are detectable, it is possible to measure a pitch 
angle. Unlike other methods, details such as distance, 
extinction, etc. do not need to be known in order to 
measure pitch angle. Additionally, the measurement of 
pitch angle is independent of cosmological assumptions. 



Upon examining t he GOODS (Great Ob servatories Ori- 
gins Deep Survey; Dickinson et al.pOOS ) fields, we have 
identi fied 224 spiral galaxies with spectroscopic (Barger 
et al | [2b08) and photometric ( |Wolf et al.||2004[ ) redshift 
[z) data for GOODS North and South, respectively. Of 
these 224 spiral galaxies, 179 galaxies lie in the range 
z < 1, 43 galaxies in the range 1 < z < 2, and two galax- 
ies with redshift greater than z = 2. So far, we have 
measured their pitch angles using the previous version 
of 2DFFT and are planning on remeasuring the sample 
using the new iterative version of the code. This work 
demonstrates that it is not uncommon to be able to mea- 
sure pitch angle for galaxies beyond a redshift of one 
Current work a lso includes artificially redshifting (Bar 



den et al. 2008) this GOODS sample of spiral galaxies 
in order to te st the completeness of the GOODS fields 
(Shields 2012). Artificial redshifting allows us to pre 



diet at what distance spiral arms are no longer visible 
and thus pitch angle is immeasurable. Results thus far 
show no indi cations of a relation ship between pitch angle 
and redshift ( Shields et al.|2010 l, but this matter will be 
further explored. 

6.4. Continuing Efforts 

It is encouraging to notice increasing interest in the 
measurement of galactic spiral arm pitch angle in the 
astronomi cal comniunity , as evidenced even in this past 
year (Kendall et al.|2011[[Mart mez-Garcfa|2012|) . Recent 
involvement from the computer science community has 
also been initiated with the use of computer visio n tech- 
niques to mea sure galactic spiral arm pitch angle (Davis 
fc Hayes|2012 ) . We feel that our modification to the pre- 
viously established Fourier methods for measuring galac- 
tic spiral arm pitch angle is a marked improvement and 
helps to satisfy a growing demand for its rapid and ac- 
curate measurement. Furthermore, we are glad to make 
the code publicly available to the benefit of anyone inter- 
ested in furthering a much-needed exploration of spiral 
galaxy structure. 
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